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I Abstract 

Data for (f) 'yiv^'^) a-re analysed using tlie KK loop model and com- 
&X). pared with parameters of ao(980) derived from Crystal Barrel data. The 

^ ' riTT mass spectrum agrees closely and the absolute normalisation lies just 

within errors. However, BES parameters for /o(980) predict a normali- 
sation for (/} 7(71'' tt") at least a factor 2 lower than is observed. This 
discrepancy may be eliminated by including constructive interference be- 
tween /o(980) and a. The magnitude required for a — > KK is consistent 



(N 



< 
(N 



X 



PACS: 13.25.Gv, 14.40. Gx, 13.40. Hq 



■ with data on tttt — > KK. A dispersion relation analysis by Biittiker, 

^ ' Descotes-Genon and Moussallam of tttt — > KK leads to a similar con- 

, elusion. Data on tttt — > rjri also require decays of a to rjTj. Four sets of 

^Nj ' TTTT KK data all require a small but definite /o(1370) signal. 

o 
cn 
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^ ■ 1 Introduction 

X 

^ ■ Data for 0(1020) li-qir^) and 7(7r'^7r'') from Novosibirsk [1-3] and the KLOE 

Qh! collaboration at Daphne [4,5] may throw light on ao(980) and /o(980), as Achasov 

^ ■ and Ivanchenko pointed out in 1989 [6]. A vigorous debate has followed the 

publication of data. An extensive list of references is given in a recent paper of 
Boglione and Pennington [7]. 

Fig. 1 shows the model conventionally assumed for radiative decays. Both 
the (j) K^K~ decay and rescattering at the /o (or ao) vertex are short-ranged. 
5^ \ In the intermediate state, kaons propagate to larger radii where they radiate a 

photon through an electric dipole moment. 

These data provide information only on the lower sides of ao(980) and 
/o(980), i.e. their coupling to yyvr and tttt. There is no direct information on 
their coupling to KK ^ except by assuming the KK loop model and using the 
absolute normalisation. For ^(rjir^), this normalisation depends on the product 
(7^[ao(980) r/Tr] x (7^[ao(980) — > KK]. In order to test the model, a comparison 
is made with parameters of ao(980) from high statistics Crystal Barrel data on 
pp — > 'K^{ri'n^) [8] and pp (jjirji:^) [9]. There is agreement within the errors, 
supporting the KK loop model. 
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Figure 1: The KK loop model for (p radiative decays. 



The BES collaboration has made an accurate determination of /o(980) pa- 
rameters from J/^* 07r+7r~ and (l)K^ K~ data [10]. Their values of pre- 
dict a normalisation for (j) — * 7(7r°7r°) a factor 2 smaller than observed in both 
Novosibirsk and Daphne data, if /o(980) alone is responsible. This discrepancy 
is resolved by adding a broad component coupling to both tttt and KK and in- 
terfering constructively with /o(980). It may be parametrised as the high mass 
tail of the a pole, although other interpretations arc possible, e.g. those of Au, 
Morgan and Pennington [11] and of Anisovich, Anisovich and Sarantsev [12]. 

This broad component should also appear in tttt — > KK. It will be shown 
that it can indeed be accomodated naturally there. A similar broad component 
is required to fit data on tttt rjr]. A combined fit is made to data on tttt elastic 
scattering, tttt KK and tttt — > rjr]. The tttt phase shifts and elasticities are 
taken from the re-analysis by Bugg, Sarantsev and Zou [13] of moments from 
Cern- Munich data [14]. 

Formulae are discussed in Section 2. Readers interested only in results can 
skip this formalism. An awkward problem arises in parametrising tttt KK 
near 1 GeV. The broad component and /o(980) both go through 90° at ~ 1 
GeV. It is difficult to find a parametrisation which accomodates both resonances 
while still satisfying unitarity. The scheme which is adopted here is driven by 
the features of the data and is explained in Section 2. 

Sections 3 and 4 then report the fits to data on </> — > 7(7777'^) and 7(7r°7r°). 
Section 5 concerns the fit to tttt KK and 7777, using parameters consistent 
with (j) 7(7r''7r°). A brief summary is given in Section 6. 

2 Formulae 

2.1 (j) radiative decays 

Close, Isgur and Kumano [15] give formulae for the mass distribution of tttt or 
ITT] pairs predicted by the KK loop model. It is important to keep track of the 
absolute normalisation and also write formulae in a way where the comparison 
between 7r°7r° and ttt] is as simple as possible. With some trivial re-arrangement 
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of symbols, their eqns. (3.6) and (3.7) may be written for (j) 7(7r°7r°) as: 



dm = IdW' 



m 



K 



^7 = • (3) 

Here s is the invariant mass squarc;d of the tttt pair, fc-y is the photon momentum, 
m are masses, a is the fine-structure constant ~ f /f 37 and /(a, 6) is a formula for 
the KK loop integral given by eqns. (3.4) and (3.5) of Close et al; the quantities 
a and h are defined there. The cubic dependence on arises from gauge 
invariance and is the familiar factor for an El transition. The denominator 
D{s) in eqn. (1) is written for /o(980) in the Flatte form: 

D{s) =M^-s- -^[gl^P^^s) + glKpKK{s)]; (4) 

the factor 1/IQ-k follows the definition used by the Particle Data Group [16] in 

their eqn. (38.17), after allowing for a factor Air from integrating the S-wave 
over angles. Eqn. (4) is the form used by KLOE. Note that values of g"^ in this 
equation refer to the sum over charge states: 

gl^ = 35^0,0 = (3/2)52^,_ (5) 

9%K = "^glcOKO = '29'k+K-- (6) 

The BES collaboration writes the Flatte form using 

»! - |j (7) 

Below the KK threshold, pkk{s) in eqn. (4) needs to be continued analytically 
as +i\J Arri^j^l s — 1. 

Eqn. (3.1) of Close et al. defines g"^ via: 

r(0 K+K-) = - ^^K?i^- (9) 

The coupling constants for /o decay are defined to have dimensions of GeV^, 
but 5^ is dimensionless and so is 7(o, 6). Dimensions of eqns. (1) and (2) balance 
between left and right-hand sides; this was not apparent in the papers of either 
Close et al or Achasov and Ivanchenko, since they replace g^j^gj^K by a single 
g'^. The two powers of g'^ appear explicitly in the KLOE publication [5]. The 
relative normalisation of data on — > 7(7r°7r°) and tt — > 7(r?7r'^) will play an 
important role: 

{dr/dm)^o^o ^ 9%kk{9%^^I'^)\D{s)\I^P^^{s) 

(dr/dm),.o gl,KK9l,.r,\D{s)\lPvM ' ^ ' 
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2.2 Identical particles 

At this point, it is necessary to pause and ask whether decays to 7r°7r° are 
affected by the fact that they are identical particles, whereas tt^t] are not. This 
is a tricky point. An J = /o couples to the isospin combination (tt+tt" + 
7r~7r~ —TT°n^)/\/3, where the coefhcients are isospin Clebsch-Gordan coefficients. 
Experimentally, the tt+tt" integrated cross section may be obtained by counting 
TT"*" over An solid angle. There are two amplitudes, one for producing a 7r+ 
at angle 0, say, and the the other for producing a 7r~ at that angle and an 
accompanying tt"*" at angle (^ + 7r). Both are in S- waves, so at angle 6 there is a 
total amplitude for tt"*" of 2/\/3. Integrating over in solid angle, the intensity of 
observed tt"*" is (4/3)47r. For the 7r°7r° case, there are likewise two 7r° at angles 
9 and {0 + tt) so the total 7r° amplitude is 27r°/\/3. In this case, the integration 
should be only over 27r solid angle, so as to avoid counting vr" pairs twice; the 
integrated intensity is (4/3)27r. The total tttt intensity is Stt, compared with the 
total r]7r intensity of 47r, where there are no effects for identical particles. 

So the TTTT amplitude is increased by a factor 2 by the identity of the two 
pions. However, it is essential to remember that this factor 2 appears also in 
the Breit-Wigner denominator. The Breit-Wigner amplitude becomes 

. ffj<:+K-g7r077o(2/V2)(l/v^) 

M^-s-i[2gl„p^^+g%KpKK]' ^ ' 

In the numerator, the factor 2 allows for the doubling of the tttt amplitude by 
identity of the two pions and the factor l/\/2 allows for the fact the tt^ are to 
be considered only over 2tt solid angle. However, experimentalists do not write 
the Breit-Wigner amplitude in this way. They ignore identical particle effects 
and write it in all cases as 

J- Gk+k-Gt^tt/ , , 

^ - M^-s-i[Gl„p^^ + G%j,pKKV ^ ' 

where G = g\f2. It is essential to remember this convention when using Breit- 
Wigner parameters deduced in the BES experiment. The upshot is that in 
comparing with Kloe data, one should use eqn. (12) and there is no effect from 
the identity of the pions in the final state. 

2.3 How to treat the a pole 

In the KLOE analysis of ^ 7(7r°7r°), a further contribution was included 
from the u pole. This requires some discussion. For elastic tttt scattering, the 
amplitude may be written 

fel - ^ (13) 

and there is a zero in N{s) at the Adler point s ~ 0.5m^, just below threshold 
[17]. The Adler zero is a basic feature of Chiral Perturbation Theory and figures 
prominently in the series of papers on a, k, /o(980) and ao(980) by Oset and 
collaborators [18-23]. Data on J/^ — > ua may be fitted in both magnitude and 
phase using 

fprod oc (14) 
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with the same denominator but a constant numerator [24]. The justification is 
that left-hand singularities due to coupling of a to wJ/^ are very distant. The 
question is whether the Adler zero is needed in the amplitude for fitting KLOE 
data or not. 

An interesting intermediate case is the production of the tttt S-wave in T' — > 
Ttttt. Recent BELLE data determine tttt mass spectra for the transitions 4S 
IS, 3S IS and 2S IS [25]. The first and last follow elastic scattering 
closely. The transition 3S — > l/S deviates from this and shows a peak towards 
the TTTT threshold, reminiscent of that observed in J/^ ^ uj<t, where the Adler 
zero is absent. The combined mass of T(3S) and 27r lies close to the BB 
threshold, and Moxhay suggests that the decay proceeds through coupling to 
BB intermediate states [26]. These data suggest that for small momentum 
transfers the Adler zero is needed in the amplitude unless there is a flavour 
change. Prom Fig. 1 for (p radiative decays, it is plausible that the right-hand 
vertex should be the same as that in the on-shell KK tttt process, since the 
kinematics are similar. It will be shown that this successfully fits the data. 

2.4 p — > 7^7 and 7r°7 

There arc also amplitudes due to </> — > 7r°p, p — > 7/7 and p — > tt°^. Formulae 
for their intensities and interferences with ao or /o are given by Bramon et al. 
[27] and Achasov and Gubin [28]; a small correction to [28] for the intensity of 
(p pn, /9 — > r/7 is made by Achasov and Kisclev [29]. Also eqn. (20) of Ref. 
[28] requires an additional factor (m| — s) in the numerator; this correction is 
made in a second paper of Achasov and Kiselev [30] in their eqn. (21). 

The pir contribution to 7(7r"7r") peaks at m^^ = 550 MeV and is small. 
Interferences with /o(980) have rather little effect on the fit to 77r°7r°. However, 
the contribution of p ^ rjTT^ is quite significant in 7(777r''). 

Branching fractions for these processes will now be discussed, starting with 
p 7r°7. Achasov et al. report r(p 77r°) = 73.5 ± 11 keV [31]. This agrees 
with the PDG average for r{p^ T""^) of 68 ± 7 keV. The weighted average 
is 70 ± 6 keV. The PDG gives an average for r((/) ^ pn + Tr+'K-Tr'^)/T total = 
0.151 ± 0.009. Dividing this by 3 for neutral final states and assuming it is all 
pTT, the branching fraction of (j) ^ pw^ tt^tt^'J is 2.4 x 10~^. 

The experimental groups apply cuts to remove (j) con^, co n^-y. It is 
hard to estimate the fraction of pn^ events which this removes; a sharp cut 
in the p line-shape at the quoted mass cut leaves ~ 40% of pw^ events, i.e. a 
branching fraction of (/) ^ 7r°7r° ~ 0.96 x 10~^. The uncertainty is not serious, 
since most events lie at low 7r°7r° masses and have little effect on the /o(980) 
peak. 

The arithmetic for the p7r° contribution to 7(r?7r'^) goes similarly. The 
weighted mean of SND [31] and CMD2 [32] values for T{p 'y'n)/Ttotai is 
(2.95 ± 0.26) X 10"^. Taking the pir^ branching fraction to be 0.151/3, 

as above, and using the KLOE branching fraction 4> ^{rj-K^) = 0.83 x 10^**, 
the ratio ((/> — > P'n'jP l''l)/{4> all "yT]n°) = 0.18 ± 0.02. These events do 
contribute significantly to the fit; interference with 700 (980) is included. 

2.5 Formulae for tttt KK 

There is a major question how to combine the amplitudes for /o(980) and the 
broad a. Data on tttt elastic scattering from the Cern-Munich experiment and 
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elsewhere show that phase shifts of /o(980) and a add to a good approximation. 
Both go through 90° near 1 GeV and the phase shift for the full amplitude 
passes rapidly through 180° just below 1 GcV. A natural and successful way of 
describing this situation is the Dalitz-Tuan prescription [33], where S-matrices 
of the two components are multiplied: 

Stotal = SaSb = VAVBe'^^^+^''^. (15) 

However, it is not obvious that it is right to multiply elasticities. 

An alternative approach is to use the K-matrix. It is conventional to add 
K-matrix elements for the two components [12]. The penalty in this approach 
is that each K-matrix pole appears where the total phase goes through 90°; the 
/o(980) is then described as a delicate interference between two poles displaced 
substantially from 1 GcV. In the present work that is inconvenient, since /o(980) 
is the focus of attention in cp decays; a formula is needed to parametrise it 
directly. Furthermore, it is not obvious that Nature chooses to add K-matrix 
elements. 

Since K ex tan 5, an alternative choice is 

Ka + Kb 
- 1-KaKb 

so that phases again add below the KK threshold. However, this alternative 
fails to fit TTTT — > KK data immediately above the KK threshold. This is 
because inelasticity of the cr amplitude grows fairly slowly above threshold, so 
eqn. (16) demands a rather elastic amplitude there. 

One cannot add inelasticities of /o(980) and a, since this leads to values out- 
side the unitary circle. The approach adopted here is pragmatic. The Dalitz- 
Tuan prescription is simple, fits the phase of elastic scattering adequately, pre- 
dicts a reasonable result for the elasticity parameter r] and successfully fits the 
data. The 77 parameter for the BES /o(980) drops precipitously to 77 = 0.27 at 
1.01 GeV and then rises again (see Fig. 11(a) below). Since /o(980) and a are 
approximately in phase, it seems likely that the effect of the a component will 
be to increase the inelasticity to a value somewhere between and the value for 
/o(980) alone. In practice, multiplying the S-matrix elements for Sn and ^22 
describes the data well up to 1.2 GeV, obtaining 5i2 from the exact relation 

2|5l2|2 = 1 + 1^331'- |^ll|'-|522|^ (17) 

where indices 1,2,3 refer to mr, KK and r]r]. Note, however, that 5*12 itself does 
not factorise into /o(980) and a components. This is a symptom of extensive 
multiple scattering. It may well be that there is no simple formula which fully 
describes strong coupling of the two amplitudes. 

Above 1.2 GeV, the appearance of /o(1370), /o(1500) and /o(1790) makes 
the treatment via the Dalitz-Tuan prescription intractable. There, amplitudes 
are small enough that unitarity corrections are less than normalisation errors, 
and it is adequate to add amplitudes. That is the procedure adopted here, 
making a smooth join at 1.17 GeV using Lagrange multipliers to ensure that 
the magnitude and phase join continuously. For the rjij channel, amplitudes are 
so small that unitarity plays no significant role, so one can add amplitudes. 

At the higher masses, /o(1370), /o(1500) and /o(1790) overlap significantly. 
The overlap makes the phase of the background below each resonance a sensitive 
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parameter. In order to achieve a satisfactory fit, it is necesssary to multiply each 
amplitude by a fitted phase factor exp(i^) and add amplitudes. 

In siimmary, multiplying rj parameters near threshold puts a lower bound on 
the intensity 92% of the unitarity limit, and reduces the normalisation uncer- 
tainty from 50 to 8%. At higher masses, the amplitude is fitted freely in terms 
of known resonance, ignoring possible multiple scattering from one resonance to 
another where they overlap. This allows maximum freedom to fit a and /o(980), 
the main object of the present work. 

An important point is that the /o(980) signal appears as a sharp dip in 
elastic scattering, because phases add. In (p radiative decays, it appears as a 
peak. In this case, it is appropriate to add amplitudes, since they are weak and 
the unitarity limit imposes no significant constraint; in clastic scattering, the 
amplitude must follow the unitarity limit up to the KK threshold. 



2.6 The 47r channel 

Both TTTT KK and tttt — > rir] cross sections drop fairly rapidly with increasing 
mass, and almost disappear above 1.9 GeV. To some extent, this is due to 
competition from the An channel. Unfortunately, there are almost no data on 
the 47r channel, so the s-dependencc of this inelasticity is poorly known. As in 
the work of Bugg, Sarantsev and Zou in 1996 [13], An phase space is parametrised 
empirically by a formula close to a Fermi function: 



= (18) 

1 - exp[A(so - s)] 

Here A = 2.85 GeV~^ and sq = 2.5 GeV^. This function approximates closely 
the onset of pp and aa final states. 



2.7 The TTTT S-wave 

The propagator of the broad <t amplitude is written in eqn. (20) as the sum 
of contributions from tttt, KK, rjrj and 47r, labelled 1-4. In the first three, the 
Adler zero is accomodated in eqn. (21) by a factor 

Ais)^{s-.SA)/{M^-SA), (19) 

where sa = (0.41 ± 0.06)m^. The factor 0.41 comes from the work of Leutwyler 
and Colangelo, using Chiral Perturbation Theory [34]. A form factor FFf(s) = 
eiip{—akf) is included into channels 2 and 3 (KK and rirj) in cqns. (22) and 
(23). Note that FFj(s)^ is multiplying the intensity of the KK channel, which 
is the quantity appearing in the width of the a. It is required to parametrise the 
rapid drop of cross sections with increasing mass. The same a is used for both 
channels; ki is the centre of mass momentum in each channel. The propagator 
is then written: 



D{s) 


= M'^ - s - gjzsub - ni{s) - i[gfpi - 


1- 9lp2 + gjps - 


^glis)] (20) 


9l 


= Biexp[-{s-M^)/B2]A{s) 




(21) 


9l 


= r.gfFFiis) 




(22) 


9% 


= nglFFlis) 




(23) 


FFiis) 


= ex.p{-aikf) 




(24) 
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9l 



z{s) = 

^sub 

m{s) 



2 + pi/rie 



1 



Pi 



1 + 

= z{s)-z{M^) 

s-A'P f MT 



{s')ds' 



is' - s)(s' - M2)' 



(25) 
(26) 
(27) 
(28) 



For other resonances, the exponential of eqn. (21) and the form factors FF are 
superfluous because the resonances are narrow; they are simply omitted, but 
the dispersive term m(s) is retained. The amphtude for tttt elastic scattering is 



Tn = 1 + 2^511 

Su = r]eM'^i5)=gfpi/D{s). 



(29) 
(30) 



In fitting elastic data, S-matrices from a, /o(980), /o(1370), /o(1500) and /o(1790) 
are multiplied together. 

The function gi{s) was introduced by Zou and Bugg [35]. It takes a very 
simple empirical form which succeeds in fitting the broad S-wave over the whole 
mass range from threshold to 1.9 GeV using just the Adler zero (taken from 
other work) and three fitted parameters Bi, and M. At the KK and r]r} 
thresholds, ra = gi-K/sln ^""1 ''s = gl-q/gl^- 

The motivation for the formalism is to make amplitudes fully analytic. This 
has been achieved below the region of strong 47r elasticity by Achasov and 
Kiselev [30] . Here the opening of the strong 47r channel is treated analytically 
for the first time. This will be important in demonstrating the presence of 
/o(1370); otherwise, the issue of analyticity turns out not to be crucial for 
present work. The function z{s) accounts analytically for the opening of tttt 
phase space pi(s); it is taken from the work of Leutwyler and Colangelo [34]. 
It has the good feature of eliminating the divergence at s = due to ^Jlfs in 
Pi(s) = \/l — Am^/s. In eqn. (27), a subtraction is made at s = M^, where 
the real part of the S-wave amplitude goes to zero. 
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Figure 2: The dispersive term m(s) for the A-k channel. 

The term m(s) is a dispersive contribution to the real part of the amplitude, 
allowing for the inelasticity in all channels. The inelasticity in 47r is large and 
gives rise to a slowly varying dispersive term illustrated in Fig. 2. It is far too 
wide to be confused with /o(1370), /o(1500) or /o(1790). The dispersive term 
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m(s) is also included into the Breit-Wigner forms for /o(1370), /o(1500) and 
/o(1790), again with a subtraction on resonance. This is important for /o(1370), 
where Ti-^ is dominant. 

One new result arises concerning the position of the a pole. Fitting it with 
the now fully analytic formulae described above moves its pole position to (506 ± 
30) - i (238 ± 30) MeV from the BBS value of (541 ± 39) - z(252 ± 42) MeV [36] . 
The shift arises from the elimination of the singularity at s = and dispersive 
corrections to the tttt channel. It is within the experimental errors, which now 
go down because of better control over the extrapolation to the pole. 
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Figure 3: The fit to KLOE data [5] where (a) rj — > 77, (b) r] — > 37r; in (a), dashed 

and dotted curves show the pn contribution and interference with ao(980); (c) 
the energy dependent factor appearing in dT/dm, but omitting |/(a, 6)p, (d) 
\I{a,b)\^ itself. 



3 Fit to liw^) 



If parameters of ao(980) are fitted directly to — > 77771° data of Ref. [5], the 
problem is that the upper side of the resonance is not visible, resulting in values 
of M, g"^^ and gj^^ which are strongly correlated. 

Most determinations quoted by the PDG have been obtained by fitting sim- 
ple Breit-Wigner line-shapes with constant width to 7777 or KK alone. These are 
of no use for present purposes, since they ignore the cusp at the KK threshold. 
A rather precise determination, not quoted by the PDG, was made in Ref. [8], 
combining data on pp ujirjir'^) and pp {riTT^)TT^ at rest. Both have high 
statistics (280K events for r^Tr^Tr^) and very low experimental backgrounds. The 
first reaction is dominated by a;ao(980). Although a simple Brcit-Wigncr of con- 
stant width was fitted, the data determine precisely the location of the ao peak 
at half-height; the error quoted on the central mass is 1.23{stat) ± 0.34{syst) 
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MeV and for the width is 0.34:{stat) ± 0.12{syst) MeV. Those constraints go 
a long way towards breaking the correlation between couplings to rjn (which 
dominates the lower side of the resonance) and KK (which dominates the up- 
per side). In titt^tt^, the ao(980) appears very conspicuously in the Dalitz plot 
as an 'edge' at the KK threshold; at this edge, the phase turns abruptly by 90° 
in the Argand diagram. The interferes with the broad a in the tttt channel. 
This interference determines rather precisely the phase variation of the Oq as 
a function of 7777 mass. On the other hand, there is some uncertainty in the 
parametrisation of the cr as a function of tt^tt^ mass. Taking into account all 
uncertainties, quoted parameters are M = 999 ± 5 MeV, = 221 ± 20 MeV 
and r = g%K/gl^ = 1.16 ±0.18. 

The fit to KLOE data with these parameters immediately reproduces the 
line-shape accurately. The fitted normalisation is low by 32%, which is close 
to the combined errors, as follows. Prom gi^^, there is an error of 19% in the 
magnitude predicted for KLOE data, and from r there is a further error of 15%. 
The error quoted by the PDG for the normalisation of ^ jrjn is 6%. To 
achieve agreement in normalisation, it is necessary to stretch all parameters by 
Icr as follows. For oq, g'^j^ has been increased to 241 MeV, and r = g^x/f^ri-rr 
to 1.34; this requires reoptimising M = 991.5 MeV. The branching fraction of 
p rjTT has been increased to 0.20 and the normalisation of KLOE data has 
been decreased by 6%. The final fit is shown in Figs 3(a) and (b). Dashed and 
dotted curves show the intensity of the p ^ rjir signal and its interference with 
ao(980). A minor detail is that a form factor exp— /3A;^ is included to allow for 
the charge radius of the kaon, 0.560 fm. [16]. 

Figs. 3(c) and (d) illustrate the strong energy-dependent factors contribut- 
ing to the intensity. Fig. 3(d) shows |/(l,6)p from the KK loop. It falls 
dramatically as the kaons go off the mass-shell. Fig. 3(c) illustrates the re- 
maining factors from mk'^pjj-^, eqn. (2). These factors inflate the lower side of 
ao(980). The agreement with the oq parameters of Rcf. [8] for both line-shape 
and normalisation confirms Achasov's KK loop model at the level of ^ 15% in 
amplitude. That conclusion will be important in considering the 77r*'7r° data. 
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Figure 4: Fit to KLOE 7(77^77°) data with /o(980) -I- pir when the normalisation 
is floated. 
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Figure 5: The magnitude of the tttt — > KK amplitude from Biittiker et al [39] 
fitted using the Roy equations and experimental data of Cohen et al [40] . 



4 Fit to ^ jiTT^TT^) 

KLOE data on ^ 7(7r°7r°) [4] can be fitted fairly well with BES parameters 
for /o(980) if the absolute normalisation is floated, see Fig. 4. The best fit is 
obtained by fitting the form factor FF freely and stretching BES parameters to 
the limits of their quoted statistical and systematic errors as follows: M = 0.957 
GcV, gl = 0.190 GcV, gj/gf = 3.75. However, = 122, compared with 60 for 
the fit shown later in Fig. 6 including a KK. The data are taken from Fig. 4 
of the KLOE publication [4] after subtracting their estimated background; they 
are not taken from their Table 5, where the background has been unfolded in 
such a way as to accomodate a dip near 520 MeV. 

There is however a serious problem with the absolute normalisation. The fit 
with /()(980) alone has a normalisation at least a factor 2 lower than the data. 
The program is identical to that used for 77771° except for trivial changes in the 
parameters of the Breit-Wigner amplitude and the change from the mass of the 
rj to the tt". A direct comparison of the relative normalisation of 777*^77° and 
7ry7r° has been made using eqn. (10). 

A variety of explanations for the normalisation are possible, but look odd 
in view of the agreement in normalisation for 7777r*' data. Firstly, Oiler remarks 
that there may be a contact term involving two if ° in the loop diagram, and a 
photon radiated from the (f) decay vertex [37]. Secondly, Oset and collaborators 
have suggested that (f) may decay through intermediate states involving K*K, 
where K* may be i4:*(890), ifi(1270) or is:i(1400) [38]. 

A further alternative is proposed here. This is that the broad a cou- 
ples to the KK loop. It will be shown in Section 5 that data on 7r7r — !■ 
KK require this broad component. Independent support is available from 
the work of Biittiker, Descotes-Genon and Moussallam [39]. They apply the 
Roy equations to fitting data on nK ttK and 7r7r KK. For the present 
discussion, this corresponds to calculating the tttt KK amplitude below 
the KK threshold as a continuation of the amplitude above threshold using 
(i) the dispersion integral over the physical tttt KK process, (ii) nearby pole 
terms due to t and w-channel exchanges of p(770) and iir*(890). Their Figure 
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12, reproduced here as Fig. 5, shows the magnitude |(7q| of their mr KK 
amplitude both above and below the KK threshold. The /o(980) signal is su- 
perimposed on a strong broad component. Achasov and Kisclcv [29] similarly 
include couphng of a to KK, but they fit parameters of /o(980) freely. 




0,4 0.6 0.8 1 
M (GeV) 




0.4 0.6 



0.4 0.6 0.8 1 

M (GeV) 




M (GeV) 



0.4 0.6 0.8 1 
M (GeV) 



Figure 6: (a) Fit to KLOE 77r*'7r° data after background subtraction; (b) con- 
tributions from /o(980) (full curve), a (dashed) and interference (dotted); (c) 
an enlarged view of the contribution from pir, p 77r°; (d) I^q] from my fit 
(full curve) and the contribution from a alone (dashed). 



For the fit presented here, the /o(980) contribution is fixed in magnitude 
to BES parameters. The charge form factor of the photon to KK is fitted 
freely. The magnitude of cr ^ KK optimises with a threshold ratio r2 = 
dKK/dvv — 0.6±0.1. Its phase with respect to /o(980) requires some comment. 
The production process is electromagnetic, and will produce a phase shift only 
of order a = 1/137. In tttt clastic scattering, both a and /o(980) have a phase 
shift of ~ 90° at 990 MeV. Their relative phase is constrained in fitting KLOE 
data within ±10°, the combined error from a phases and the mass of /o(980). It 
fits naturally at the edge of this band; if the constraint is removed, changes 
by only 2 and the fit hardly changes. 

The resulting fit is shown in Fig. 6(a). Panel (b) displays the /o(980) 
intensity as the full curve, the a as the dashed curve and the interference between 
them as the dotted curve. Panel (c) shows the fitted pn signal (full curve) 
compared with data at low masses. Its magnitude is fixed according to the 
arithmetic of Section 2.4. A marginal improvement may be obtained by fitting 
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it freely, but there is no significant effect on the fitted a amphtude. 

Fig. 6(d) shows \g%\ from my fit to tttt — > KK as the full curve. The peak 
due to /o(980) is slightly lower than that in Fig. 5 because the /o(980) of 
BES is slightly broader than that of Biittiker et al. Secondly, the fit requires a 
fairly strong form factor to cut off the lower sides of both /o(980) and a. It is 
parametrised as cxp — 8.4A;^,, where k is in GeV. Comments on the origin of this 
factor will be given in Section 4.1. 




0.8 1 0.5 0.75 



Figure 7: Fits to Novosibirsk data on ^ 7(r/7r°) [2]''kifif'$(7r°7r°) [3]. 

For completeness, Fig. 7 shows fits to Novosibirsk data on ^ — > 'y{ri'iT^) and 
77r''7r°. The point in 7(r]7r") at 0.98 GcV is clearly too high and is slightly 
inflating the PDG average for the branching ratio </> 7(7771"). 

4.1 Sensitivity to fitting parameters 

The heights of the peaks at 970 MeV in both ^{itk°) and 7(7r'^7r°) have lim- 
ited sensitivity to resonance parameters. In the denominators D{s), g^^p^^ or 
drt-nPn-^ appear in the imaginary part, while M and g\x \J 4Af^ /s — 1 appears 
in the real part. There is considerable compensation between the numerator 
N[s) and denominator D(s) of the Breit-Wigner amplitudes. For 7(r77r"), it is 
necessary to constrain the location of the ag peak at half-height within, say, 2 
or 2.5 standard deviations of the values quoted in Section 3, adding statistical 
and systematic errors linearly. With this constraint, M, g^^ and g\x ma-y be 
varied freely. Because of the feedback between iV(s) and £'(s), it is only just 
possible to accomodate the 32% discrepancy by including also flexibility in the 
pTT contribution and uncertainty in the KLOE normalisation. 

For /o(980), the BES tttt peak must again be constrained in mass and the 
full width must be constrained to 34 ± 4 MeV. There is more flexibility in 
9kk-i since the BES /o(980) KK signal interferes with a background due to 
cr KK; the error on g\x-> namely ±0.25(siai) ± 0.21(syst) quoted by BES, 
represesents one standard deviation. It turns out that with this flexibility, the 
absolute normalisation of the /o(980) peak cannot change by more than ±15%. 
Some fits reported in the literature have changed g\^ and g\xl An factors 
of 2; that is unrealistic. 

The peak at 970 MeV in KLOE data determines the magnitude of the back- 
ground amplitude there. The fit is not unduly sensitive to the parametrisation 
of TTTT phase shifts. What matters for 7(77° 7r°) is only the way the a and /o(980) 
amplitudes go out of phase ofi^ resonance. Any form which reproduces tttt phase 
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shifts within errors of, say, ±4° is adequate. The fuU apparatus of eqns. (18)- 
(30) has been used to scrutinise the combined fit with tttt — > tttt, KK and 777/. 
An equivalent and more convenient form which fits all the data equally well 
may be obtained by (i) dropping the dispersive contribution m(s) completely, 
(ii) also the tiny contribution from the 47r channel below the KK threshold, (iii) 
replacing the effect of m(s) for KK and 7777 below their thresholds by a form 
factor exp(-5.2|A;|.|), with kx in GeV/c. Then M = 0.7366 GeV, = 0.6470 
GeV^, B2 = 4.271 GeV^, B4 = 0.00221 GeV^. Values of other parameters 
will be discussed below but will be collected here for convenience of reference: 
ra = 0.20, = 0.6. 

A crucial parameter in fitting a — > KK in KLOE data is the form factor of 
Eqn. (24), FF = exp(— 8.4fc^), where is in GeV. Without this factor, the 
background amplitude is too large at low mass; one can see this in the difference 
between Figs. 5 and 6(d). This form factor has at least three possible origins. 

Firstly, it may be due to the size of the KK cloud. This would require an 
RMS radius of 1.4 fm. It is well known that a source of finite size produces a 

form factor sin(fcr)/fcr = 1 — k^r^/6 + The second sheet pole of /o(980) lies 

at (998 ±4) - 7(17±4) MeV. The binding energy is so close to the KK threshold 
that the /o(980) inevitably has a KK cloud of roughly this size, resembling the 
long-range tail of the deuteron. The pole for ao(980) is much further away, at 
1032 — 785 MeV, so the RMS radius of ao(980) has conventional dimensions and 
a weaker form factor. 

A second straightforward possibility is that /o(980) and the broad a com- 
ponent mix over the mass range where they are both large. This mixing will be 
confined naturally to their region of strong overlap. 

A third possibility is that |5q| fitted by Biittiker et al. has some flexibility 
in the range s = 0.2 — 0.8 fm. It is well constrained near 1 GeV, and also for 
s = —0.5 to GeV^ (by data on nK ttK). In between, Adler zeros in tttt at 
s ~ 0.5to^ and in coupUng to KK at s ^ O.bm'j^ may pull j^ol down. Data on 
(j) 7(7r°7r°) appear to be the only data directly sensitive to this possibility. 

5 Fits to TTTT KK and rjr] 
5.1 Introduction 

A problem with data on tttt KK is that cross sections from different exper- 
iments differ in normalisation by up to a factor 2. This problem is minimised 
here by using the BES parameters for /o(980) to assist the absolute normali- 
sation. All TTTT — !■ KK data are rescaled to the BES normalisation. However, 
one can still not rely on tttt KK data to separate branching fractions of 
/o(1370), /o(1500) and /o(1790) to KK and rjrj. Vastly better determinations 
are available from the analysis of Crystal Barrel data on pp annihilation at rest 
to 377°, ?77r°7r°, tji^k^ and several KKtt charge combinations. All of these data 
have been fitted self-consistently by Sarantsev and collaborators; decay widths 
to TTTT and rji] arc taken from this analysis [41]. Their errors are used to set lim- 
its within which parameters may vary; Table 1 shows these limits and values at 
which they optimise. Note however, that the present data really only determine 
the products Ftttt x ^kk and F^jr x F^^. Values of F determined from 

F47r = ^totai — Etttt — ^KK — ^rjri at each resonancc mass. 

Some comments are needed on /o(1370) and /o(1790). There has been some 
unnecessary controversy recently concerning the existence of /o(1370), so its 
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State 


/o(1370) 


/o(1500) 


/o(1790) 


M 


1290-1335 (1290) 


1495-1510 (1501) 


1770-1800 (1800) 


^ total 


200-280 (280) 


110-135 (135) 


180-280 (280) 


V 


34-58 (48) 


35-39 (35) 


30-100 (72) 


'^Kk/^tztz 


0.04-0.235 (0.235) 


0.22-0.272 (0.272) 


0.16-1.2 (0.47) 




0.08-0.18 (0.18) 


0.04-0.18 (0.11) 


0.1-1.0 (0.20) 



Table 1: Limits used in the fit (in MeV); optimised values are in parentheses. 



properties will be outlined briefly; it will be the subject of a separate publication 
giving full details from latest analyses of Crystal Barrel data. It appears most 
clearly in pp — > Btt*^, where it is statistically at least a 50 standard deviation 
signal. Statistics are 600K events, with almost no experimental background. In 
early analyses of these data, there was some sensitivity to the parametrisation of 
the TTTT S-wave amplitude, with which it interferes. The behaviour of the broad 
S-wave amplitude is now known much better. 

The /o(1370) decays dominantly to An with r(7r7r)/r(47r) < 0.2. The An in- 
elasticity increases rapidly with mass, so it is essential to treat the s-dependence 
of this width fully. The result is that the tttt signal peaks at 1320 ± 30 McV, as 
in the analysis of Anisovich et al. [41]. The 47r channel peaks approximately 75 
MeV higher, because of the rapidly expanding phase space. The apparent width 
in TTTT and 47r channels is likewise different. PDG listings assign unreasonably 
large errors because these facts have not been taken into account in most anal- 
yses. The s-dependence of the 47r width is taken into account here, including 
the associated dispersive correction to m{s) in the Breit-Wigner amplitude. 

In elastic scattering, the intensity of /o(1370) is < 4% of the unitarity limit. 
It is swamped by the /2(1270), which has a spin multiplicity 5 times larger. 
Achasov and Shcstakov point out a rather large error in for /2(1270) [42]. 
This adds to the difficulty of normalising Cern-Munich moments in the region 
of /o(1370). 

The /o(1790) is required as well as /o(1710) by recent BES data. The 
/o(1710) appears clearly in J/^P coK'^K" [43]. It is conspicuously absent 
in J/^ WTT+TT" [36]. Those two sets of data require tttt/KK < 0.11 for 
/o(1710) with 95% confidence. In J/^ (pTT^TT~ , there is a distinct rnr signal 
at 1790 MeV, but no significant signal in (pK~^K~ [10]. The branching frac- 
tion TTw/KK is a factor 25 larger than allowed for /o(1710), hence requiring a 
separate resonance /o(1790). There was previous evidence for it in J/^ j4tt 
[44,45] and pp rjrjn'^ [46]. Here, both /o(1710) (with its upper limit for tttt 
coupling) and /o(1790) (with its upper limit for KK coupling) are included. 
However, the present data do not distinguish cleanly between them. 

5.2 TTTT — r]r] 

Fig. 8 shows fits to tttt KK data and rfrj. The fit to rjrj in panel (d) will be 
discussed first, since it is simple and concerns only one set of data, from GAMS 
[47]. The normalisation is such that the unitarity limit is at an intensity of 
0.25; this arises from the relation T12 = 6*12/22 between T and S-matrices. The 
full curve shows the overall fit, including contributions from a (dotted curve), 
/o(1370), /o(1500) and /o(1790)//o(1710). The main features are a gradually 
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falling a ^ rjr] signal and destructive interference with /o(1500), which is mainly 
responsible for the dip at 1.45 GeV. The peak at higher masses comes from 
/o(1500) interfering with /o(1710) and /o(1790). The Argand diagram is shown 
in Fig. 9(b). 




M(GeV) M(GeV) 

Figure 8: (a) Fit to tttt KK (full curve); circles show data of Lindenbaum 
and Longacre [48] and crosses those of Martin and Ozmutlu [52]; the dashed 
curve shows the intensity of /o(980) and the dotted curve that from a; (b) fit 
to phases of Etkin ct al [54]; (c) fit to the data of Cohen et al. [40] (squares) 
and Polychronakos et al. [53] (triangles); (d) fit to data of Binon et al [47] for 
TTTT — > r}r] (full curve); the dotted curve shows the a contribution. 

The fitted ratio = g^/gf at the ryry threshold is 0.20±0.04. The error arises 
from normalisation uncertainty and a possible weak coupling of /o(980) to ryry; 
an upper limit on this coupling is obtained from BES data on J/tp ^ (/)7r+7r~ 
and <pK~^K~, from the fact that no drop is observed at the r]T] threshold in the 
TTTT mass spectrum. This upper limit is ra ~ 0.33. 

5.3 TTTT KK 

It is first necessary to review the many available sets of data. It is important 
to realise that an analysis of moments was required in order to separate spins 
J = 0, 1 and 2 and small amounts of higher J. The papers make clear that 
there is some cross-talk with the high mass tail of p(770) and with /2(1270), 
whose branching ratio to KK is not particularly well known. 
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Figure 9: Argand diagrams for (a) tttt KK and (b) tttt — > 7777; masses are 
marlied in GeV. 



Costa ct al. [50] produced data on 7r~p K^n, but found 8 alternative 

solutions, because of ambiguities over the exchange process and the possibility 
of both 7=1 and contributions to the final state. Cason et al. [51] produced 
data on tt~p K^K^ and observed a shoulder at 1300 MeV. They obtained 
two solutions, one with a narrow S-wave resonance (width ^ 95 MeV) at this 
mass and the second with a more slowly varying amplitude similar to today's 
/o(1370). The situation clarified further with data from Pawlicki et al. [49] on 
both iT~p K^K^n and 7r+n K''K^p. These data favoured the second 
solution of Cason et al. and assigned / = 0. A meticulous analysis of the t- 
dependence of data was made by Martin and Ozmutlu [52] . This analysis showed 
that TT exchange dominates at low t, as expected. The paper of Polychronakos 
et al. [53] on tt~p nKgKg is a full length paper on the data of Cason et 
al. The publication of Cohen et al. [40] gives a revised partial wave analysis 
of the data of Pawlicki et al. and again favours a structure close to today's 
/o(1370). Etkin et al. [54] reported data on ■jT~p — > KgKgn. Further statistics 
were added in the paper of Lindenbaum and Longacre [48] . 

It is noteworthy that all the later sets of data [48], [40] and [53] observed a 
definite small bump at 1300 MeV, which was christened the e at the time and 
is now /o(1370). All also observed a threshold peak due to /o(980). However, 
there are discrepancies concerning the height of the /o(980) peak with respect to 
the 1300 MeV mass range. The analyses of Martin and Ozmutlu and Longacre 
et al. agree well. The revised analysis by Cohen et al. of the Pawlicki et al. 
data gives a distinctly lower /o(980) than that of Longacre et al. And the data 
of Polychronakos et al. give an even lower /o(980) peak. 

To clarify the situation including today's information about /o(1370), the 
four sets of data from Lindenbaum and Longacre, Cohen et al., Martin and 
Ozmutlu and Polychronakos et al. were included into the inital stages of the 
analysis reported here. This was despite the fact that those of Martin and 
Ozmutlu and Cohen et al. are based on the same actual data before the moments 
analysis. The objective is to see whether the analysis favours one or the other. 
This is not the case and the fit goes midway between the two, so they are both 
retained in the final fit, but with a weight half that of Longacre et al. 

The final fit is shown in Figs. 8(a) and (b). The /o(980) produces the 
threshold peak, but its tail at high mass cannot account for the remaining 
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features. Some a KK is definitely required and the dotted curve shows the 
optimum fit. It requires a form factor FF = exp(— 5.2fc^) with k is GeV/c. 
However, some of this form factor may reflect the effect of competition between 
KK and 47r channels. The rjr] threshold contributes a drop in the intensity of 
0.02 between 1.09 and 1.18 GeV. The /o(1370) helps fit the bump at 1300 MeV. 
Destructive interference with /o(1500) is required to fit the dip at 1550 McV. 
The fit to the 1550-1900 MeV intensity is not perfect because of difficulties in 
reproducing the large phase variation reported by Etkin et al. [54]. The fit to 
their phases is shown in panel (b). However, one should be aware that other 
phase determinations show systematic differences with their data: the whole 
curve can move up or down by up to 30°, although the trend of the phase with 
mass remains similar to Fig. 8(b). These differerences probably arise from the 
way spin 2 is fitted in different analyses. Etkin et al. included /2(1525), but 
at that time the existence of /2(1565) was not known, and could lead to small 
systematic shifts. Another possible source of systematic error is that the present 
fit ends at 1900 MeV; there may be some contribution from the known /o(2020) 
which is presently not fitted. 



>, O.OB 
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Figure 10: Fits without /o(1370) for (a) tttt, (b) 777?. 

The optimum fit gives r2 = O.GIqj, eqn. (22), in close agreement with the 
KLOE data. This value can depend systematically on the s-dependence fitted 
to /o(980) and on An inelasticity. The /o(980) is fitted with a conservative form 
factor exp(— 2.7fc^), where k is kaon centre of mass momentum; this corresponds 
to a radius of 0.8 fm for the tttt KK interaction, i.e. a conventional radius. 
Fig. 8(c) compares the fit with intensities derived by Cohen ct al [40] (squares) 
and those of Polychronakos et al. [53] (triangles). The latter are systematically 
low below 1.15 GeV, so they are omitted from the final fit. 

Fig. 10(a) shows the fit without /o(1370). It fails to fit the dip below 1300 
MeV. Since that structure is observed in 3 experiments and 4 sets of data, 
it is hard to deny the presence of /o(1370); its fitted mass and width agree 
within errors with the best determinations from Crystal Barrel [41] and BES 
[10]. There is presently no alternative explanation of the 1300 MeV bump. The 
onset of 477 inelasticity is far too slow to account for structure with the width 
observed around 1300 MeV, and is anyway taken into account fully in the present 
analysis. Fig. 10(b) also shows the change in fitting tttt rjrj when /o(1370) 
is removed; however, the points in Fig. 8(d) fluctuate around the fit by more 
than statistics, so it is difficult to estimate the significance of /o(1370) — > 777/. 
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Figure 11: (a) Contributions to elasticity r] for tttt elastic scattering from /o(980) 
alone (full curve), a (dashed) and their product (dotted), (b) Fitted tttt phase 
shifts; points show data of Pislak et al. [55] below 400 MeV and Cern-Munich 
data above 600 MeV [14]. 



Fig. 11(a) illustrates the way elasticities of /o(980) and a combine. The 
dotted curve shows the result obtained by multiplying /o(980) and a contribu- 
tions, according to the model assumed here. Fig. 11(b) shows the phase shift 
fitted to TTTT elastic scattering; errors above the /o(980) are typically ±15°. 

6 Summary 

Data on — > 7(r/7r) agree within errors in both absolute normalisation and mass 
spectrum with the parameters of ao(980) from Ref. [8]; the KLOE data suggest 
a normalisation which can be accomodated by shifting ao(980) parameter to 
M = 991.5±4 MeV, = 0.241±14 GeV^, q^k/oI^ = 1.34±0.13. The errors 
quoted here are reduced compared with Ref. [8] after the combined fit with 
KLOE data. The agreement with KLOE data suggests that the KK loop model 
is reliable. If BES parameters for /o(980) are used, data on ^ — » 7(7r"7r") have 
a normalisation higher than the fit by at least a factor 2. This discrepancy may 
be solved by allowing constructive interference between /o(980) and a KK. 
The fit requires a ratio r2 = g^{a KK)/g'^{a tttt) = 0.6 ±0.1 at the 
KK threshold. Data on tttt KK also require coupling of a to KK with the 
same value of r2 within experimental errors. Data on tttt rjrj require a ratio 
r3 = g'^{a r]r])/g'^{a rjri) = 0.20 ± 0.04 at the ryry threshold. 

The interpretation of branching ratios in terms of models for /o (980) , ao (980) , 
(7 and K is discussed in an accompanying publication. 
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